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ABSTRACT 

o, 

w ■ We have developed a new inethod(E-HOLICs) of estimating gravitational 

"^ ■ shear by adopting an elliptical weight function to measure background galaxy 

' ' images in our previous paper. Following to the previous paper where isotropic 

^ ■ Point Spread Function(PSF) correction is calculated, we consider an anisotropic 

PSF correction in this paper in order to apply E-HOLICs for real data. A an 
example, E-HOLICs is applied to Subaru data of massive and compact galaxy 
clusters A370, and is able to detect double peaks in the central region of the 
cluster consistent with the analysis of strong lensing. We also study the sys- 
^ ' tematic error in E-HOLICs using STEP II simulation. In particular we consider 

i/^ • the dependences of signal to noise ratio "S/N" of background galaxies in the 

[-» ■ shear estimation. Although E-HOLICs does improve systematic error due to the 

■ ellipticity dependence as shown in part 1, a systematic error due to the S/N 

dependence remains, namely E-HOLICs underestimates shear when background 
galaxies with low S/N objects are used. We discuss a possible improvement of 
^ ■ the S/N dependence. 

1. Introduction 

Weak gravitational lensing has been widely recognized as a unique and very powerful 
method to study not only mass distribution of the universe but also the cosmological param- 
eters (see for example, Mellier 1999, van Waerbeke & Mellier 2003, Schneider 2006, Munshi 
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et al. 2008). In particular, cosmic shear has attracted much attention recently because of its 
potential to determine the so-called cosmic equation of state, namely the relation between 
the effective energy density and pressure of the dark energy. There are some detections of 
cosmic shear(Bacon et al 2000; Kaiser et al 2000; Maoli et al 2001; van Waerbeke et al 
2001: Refregier ey al 2002; Bacon et al 2003; hamana et al 2003, Casertano et al 2003; van 
Waerbeke et al 2005; Massey et al 2005; Hoekstra et al 2006). However the lensing signal of 
cosmic shear is very small and thus highly accurate shear measurement is required. The most 
popular method of shear estimation is given by Kaiser et al. 1995(called the KSB method: 
see also Luppino & Kaiser 1997; Hoekstra et al. 1998; Viola et al. 2011) where change in the 
moments of galactic light distribution by lensing is extracted from the measurement. Other 
methods of shear estimation have been also developed (Bernstein & Jarvis 2002; Refregier 
2003; Kuijken et al. 2006; Miller et al. 2007; Kitching et al. 2008; Melchior 2011). In 
relation with our E-HOLICs method, we mention the work by Bernstein & Jarvis 2002 and 
Melchior 2011 who also introduced an elliptical weight for the shape measurement. All of 
them are very challenging and successful for cluster lensing. Although some of them have 
already been used for cosmic shear, it is argued that none of them have achieved the required 
level of accuracy to measure the cosmic equation of state in a few percent level. 

On the other hand there are planned cosmic shear observations in very near future for 
the purpose of measuring the cosmic equation of state. In this situation, it is urgent to 
develop sufficiently accurate method of shear estimation. In our previous paper(Okura & 
Futamase 2011, we call part 1) we have developed a new method of the evaluating shear 
which is based on KSB method by adopting an elliptical weight function to measure the 
shape of background galaxies. Our method is a natural development of our previous studies 
of weak lensing analysis which uses higher order multiple moments of shape of background 
galaxies. We called the method as HOLICs(Higher Order Lensing Image Characteristics) 
which are quantitates with a definite spin properties made of higher-order multipole mo- 
ments(Okura et al. 2007). We have shown that Oct-HOLICs is an unbiased measure of 
Flexion and is successfully applied to a galaxy cluster Abell 1689 to reveal substructure 
in the central part of the cluster(Okura et al. 2008) and spin-2 HOLICs which is spin-2 
combination of higher order moments increases the lensing information of image and is a 
useful method to apply cosmic shear measurement (Okura & Futamase 2009). We called 
our method as generalized as elliptical weighted HOLICs(E-HOLICs) method because of the 
elliptical weight. We have showed that HOLICs method with an elliptical weight is able to 
measures the lensing distortion more accurately by weighting highly to brighter region of 
image than in the standard KSB method, and thus it can reduce effects of systematic error 
and random noise more effectively than KSB method. We have also calculated isotropic PSF 
correction in part 1 and tested its accuracy using STEP 2 simulation data. Following to part 
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1, we present a method of correcting PSF containing anisotropic part in this paper. This is 
necessary to apply the method to real data. In the case of isotropic PSF correction, it was 
shown that the corrected shear can be obtained by solving a polynomial equation. However, 
since anisotropic part has an arbitrary direction, it is natural to expect that the shear after 
PSF correction is obtained by solving coupled two polynomial equations. We succeed in 
reducing these coupled equations into tractable forms by dividing anisotropic part of PSF 
into parallel part and orthogonal part against the direction of shear (see section 2 in detail). 
We test E-HOLICs method with STEP2 simulation data. In part 1, we can only test the 
data with isotropic PSF, but in this paper we test all PSF set and obtain more detailed 
results. After establishing PSF correction, we are able to apply our method to real data of 
galaxy cluster A3 70. 

The organization of the paper is as follows. First we give a brief review of weak lensing 
and E-HOLICs method in section 2. Then in section 3, we describe the correction of PSF with 
anisotropic part in E-HOLICs method by dividing anisotropic part of PSF into anisotropic 
part into parallel part and orthogonal part. In section 4 we tested E-HOLICs method by 
using STEP2 simulation data, in particular we investigate the dependences of S/N and size 
of background galaxies in the shear estimation. Then we apply E-HOLICs method to real 
data of galaxy cluster A 370. It will be shown that E-HOLICs is able to detect two peaks in 
the central part of the cluster consistent with strong lensing analysis. Finally we summarize 
our results and give some discussion in section 6. 



2. Bases of Weak Lensing and Elliptical weighted HOLICs 

In this section, we present bases and definitions for E-HOLICs method. More details 
can be found in part 1. 



2.1. Notations and Definitions 

Here, we describe briefly our notation and the concept of Zero image. The relation 
between zero image, source and observed image will be shown in FigHJ The bases of Weak 
Lensing can be found, for example, in Bartellman and Schneider 2001. 

We use complex notation for angular positions (e.g. 9 = 9i + 162 in the image plane and 

[3 = I3i + i(32 in the source plane). 

For notational simplicity, we set the centroid of image as the origin in our coordinates 
(therefore centroid position is ^ = ^1 + 2^2 = and (3 = f3i + 1^2 = 0). We introduce an 
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imaginary plane called zero plane where shapes of all sources are perfect circles, and regard 
the intrinsic shear as the result of an imaginary lensing distortion. The sources in the zero 
plane are called zero images. This plane is introduced in order to define naturally elliptical 
window in measuring shapes of background galaxy images. Then the lens equation gives the 
following relation between the displacement in zero plane (/3) and image plane(^) as shown 
in part 1. 

~p = [i-K){e-g^e*), (1) 

where k is a dimensionless surface mass density and g^ = g^' + igr^ is the " combined shear" 
defined as 

where g^ is the shear induced by lensing and g^ is the intrinsic reduced shear. Using the 
combined shear, we can divide lensing shear analysis into two steps as measuring g*-^ from 
each background galaxies and then determining lensing shear by statistical averaging such 
that 

Henceforth, because a purpose of E-HOLICs is measuring g*" , we notate g = g'" for 
simplicity. The complex distortion is defined as 

^ = ~r> — 2' 

and the absolute values of them are notated as 



(4) 



9 = \9\ (5) 

S = \d\. (6) 

We define complex moments of brightness distribution I{6) measured with weight func- 
tion which has ellipticity S as 



a^ 



where 



»£ = 9^ 



,y I'^OjliM^'j ^ j-(ll2-Ro[«-»3))/,= (g) 
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and 0" is arbitrary scale (usually typical scale of each object). Then HOLICs is defined as 



nUi,zk,d)^^M^- (10) 



3. PSF correction in E-HOLICs method 

In this section, we give a formulation for the correction of PSF with anisotropy in E- 
HOLICs method. Generally, PSF will be arbitrary function, but we assume PSF has only 
elliptical part. 

First we define /^ as a brightness distribution of a lensed image which has ellipticity S 
and /^ is a brightness distribution of a zero image, so 

1^0) ^ I\9) (11) 

Because zero image is defined as not having ellipticity with circular weight function, we can 
obtain 

Hl{l\Zl6) = 6 (12) 

?^2(/^,Z2,0)=0, (13) 

where let 5" is a scale of weight function in the zero plane and a is slightly different from 
that in the image plane a, because lensing changes scale of metric. Next, we consider J*^° 
which is smeared I^ by circular PSF P*''°(/3), therefore we have the following relation 

r"{P) = I d^i^{^p)P'"'{(3 - ^p), (14) 

where let P{0) as a PSF function of P*'^°(/3) in the image plane, so 

P{e) = P''"0) (15) 

and because I^ and P*''" don't have ellipticity and this relation is same as eq. (ll2p and eq. 013p . 
we can obtain 

nUP,Zl5) = 5 (16) 

nl{p'"',zlo) = 0. (17) 

And we define I^ as a lensed image from P'^°{/3), so 

I^{9) = r°{f3). (18) 
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Because /**°(/3) is made by convolution of circular functions, /*'"'(/3) also doesn't have ellip- 
ticity, so we obtain 

nl{i\zl5) = 6 (19) 

Hl{r^\ZlQ) = Q. (20) 

Finally by writing eq.( !T^ in the image plane with eg. ( ITT]) . eq. (TT5|) and eq. (TT8|) . we obtain 

F{9)= f d^ipi^{iP)P{e-iP). (21) 



These equations mean that if the image (e.g. /^) and the smearing function (e.g. P) have 
same ellipticity, the smeared image (e.g. /^) also has the same ellipticity with them. Then 
the S can be determined by observing E-HOLICs, eq. fl2Tl) using /^, eq. flTQj) . On the other 
hand the observed image is the smeared /^ by real PSF P{9), 

r^'{d) = f d^ijI^{ij)P{9 - 7/-). (22) 

Thus by finding a transformation from the observed real PSF P{6) to P{6), S can be deter- 
mined. This is what we are going to do in the below. 

In KSB method, PSF is divided into anisotropic part and isotropic part. However, 
E-HOLICs method divide PSF into orthogonal part and parallel part against the direction 
of lensing distortion. Since the direction of PSF ellipticity is not always the same with the 
direction of lensing distortion, it is natural to expect that the equations of correcting PSF are 
two dimensional form. By dividing PSF into orthogonal and parallel parts, PSF correction 
can be written in two simple one dimensional forms as shown below. 



3.1. General form of PSF correction 

Here we demonstrate general form of anisotropic PSF correction. 

Let q{9) be a part of anisotropic part (it is not necessary to be whole of anisotropic 
part), P{9) can be divided into as 

Pi9)^ J£i;P,i9-i;)q{i;). (23) 

This decomposition is one of the basic assumption of KSB method, but there have been some 
discussion on the validity of this decomposition (Kuijken 1999). More careful study may be 
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necessary on this problem. Therefore /°*^ is 

r'^(^) = [ d^^l\^)P{9 - ^) = / d'i^I^ij) f d^P,{e -i,- 0)q(0) 

= Jd'<pi,i9-<P)q{<P), (24) 

where 

i,{e)^ j d^i\i,)p,{e-i,). (25) 

Let define /(^) = ^mW^((^o "~ ^-^ f'^*^!]) Z*^^)) ^'^'^ by calculating moments, we obtain 

d'^ef{d)F^'{e) = [ d'^9f{9) f d^iPI^{^)P{9-iP) 

d^ef{e) j d^i,{e-ci>)q{(^) 

rfV /"rfV/,M/(v^ + 0)g(0). (26) 



By expanding /((/? + 0), moments of 1°^^ can be expressed by combination of moments of Iq 
and q. Especially, Z| and Z^ are obtained as 



Zi{F'\d) ^ QZ^il^S) 



+Qql 



+Qql* 



ZO - -X (4^2 - 56* Zl) + -i^ {Z^ - 16* Zl + (J^^Z^^) 



^° 2^2 
1 



2^4 



:/.,'5) 



1 



2^2 



(-5Z,^) + — {Zl - 25Z,^ + 5^Zo^) 



2a4 



(/., 5) 



(27) 



Zo^(r^^5) ^ gz2(/,,5) 



where 



+^^^2 



+^^2* 



•^ (2Zr - 35*Zo^) + ^ (Z^^* - 25*Zo^ + <5*^Z,^) 



1 



1 



2a2 



(2Z,^ - ?,6Zl) + — (Z,^ - 2,5^0^ + 6'Zt) 



2(t4 



(/., ^5) 
a.,'^), (28) 



9^;^/ 






Q = fd^eq{9)wf^ 



(29) 
(30) 



andgO = l,q^ = 0(iV^O). 



3.2. Parallel and orthogonal decomposition of PSF 

Here, we introduce the concept of paraUel part and orthogonal decomposition of PSF. 
This allows us to calculate PSF correction with two simple forms. The relations in PSF 
correction is shown Figj2l In KSB method, PSF is divided into isotropic PSF P**°(^) and 
anisotropic part q°-'^^^°{0'^ as follows. 

P{d) = I d^i)P''°{iP)q''''''°{e - ij), (31) 

Instead of the above decomposition, we decompose PSF into parallel and orthogonal part 

as, 

Pi9)^ J£i,P//^ij)q^i9-^). (32) 

where the ellipticity of P//{0) is parallel to the ellipticity of the lensed image, and the ellipticity 
of qy;{0) is orthogonal to ellipticity of the lensed image. The moments of these parts are 
defined as 

q>Si ^ q; — ^^ (33) 

Qx ^ j SeqJ,e)W (^^^ (34) 

This decomposition might have the same problem with KSB decomposition eq. (l23|) . It 
remains to be seen if the higher order correction may correct the real PSF in a reasonably 
accurate level. 

For the moment we assume the above decomposition and define "parallel image" I //id) 
which is seared by P/fO) as 

1,^9)^ j<mi\ct^)P,ie-ct^). (35) 

We define qs{d) which relates P/fd) with P as follows 

P{e) ^ J £<pP//i<P)qsi9 - (j)) (36) 

and moments are defined as 

Jd^eeUsio)wm 

QsM ^ 7^ ^^ (37) 



Qs ^ I d'9qs{9)W (^^y (38) 
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where qs\ have a same direction of 6. The anisotropic part qr'^"*"" in KSB decomposition has 
an elhpticity which has the same direction of 6. 

A catalog of these definitions can be seen in Appendix \K\ 

For notational simphcity we set the direction of distortion as real, namely 

8 = 5 (39) 

Then the following expressions are obtained by definitions. 



9x2 



Qsl- 



(40) 
(41) 



3.3. Orthogonal PSF correction(determining a direction of distortion) 

Here, we present a method of orthogonal PSF correction. After correcting orthogonal 
part, I// has the direction of distortion, therefore this is same as determining a direction of 
distortion. 

From eq. fl27|) . eq. fl28l) . eq.f l33|) and eq. fj40|) . the correction of orthogonal part of PSF is 
written as 



zi{r'^,6) ^ Q.ZiiI//,5) 



+iQxq: 



^ (4Zo^ - 56* Zl) + ^ {Z^ - 26* Zl + 6*'Zt) 



(I//, S) 



+ {iQxqxiy 
Q>cZl{I//,6) 
+iZl{Ii/,6)Qyq, 



^A-^Zl) + ±,{Zt-26Zt + 6^Zt) 



{I/h S) 



H'o - ^ (^0^ - SHI) + ^ (1 + 6^) {H^, - Ht) 



2a4 



{I//,Zl6) 



= QXi^h S) [Hl{I/h Zl 6) + 2P>f (///, 5)gx2 
Zl{r'^,6) ^ QxZl{I/,,6) 



(42) 



~r^Vx'?X2 

+ («Qxgx2)* 
QxZl{I//,6) 



2^2 



{2Zl* - ?>6*Zl) + —^ {Zt - 26* Zl + 6*'^Z\) 



2a^ 



2^2 



{2Zl - UZl) + -ij [Z\ - 26Zt + 6'Zt*) 



2a^ 



(///, S) 
{I/h S) 



(43) 
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where we used the condition that all moments of right hand have only real parts. Therefore 
we obtain 

Hl{r'',5) ^ Hl{I„,Zl5)+iP^{I,/,5)q^^, (44) 

Thus we can divide the observed ellipticity into real part and imaginary part as follows. 

Re[nl{r'\5)\ = HliI//,Zl6) (45) 

lm[nlir'',6)] = P^^il//,6)q.l. (46) 

Similarly the moments of star J°^** can be divided into as 

Re['Hlir"''',6)] = Hi{I^,Zl6) (47) 

lm[Hl{r'^''^,6)] = P^{I^,S)q^l (48) 

Thus we have 

'^^" Px^(///,5) " Px^(/;,5) ■ ^ ^ 

In real analysis, first we must find a direction which satisfies the second equality of eq. (H9l) . 
and the removal of imaginary part in eq.(jH]) corresponds to the correction of orthogonal 
PSF. 

Let's consider that the observed ellipticity of image is A + iB and star is a + z6 in 
arbitrary basis, respectively. The direction of distortion (ps is determined by the following 
equation. 

_J p,^(r,,s)B-p^(i,,s)h \ 



a 



3.4. Parallel PSF correction(determining absolute value of distortion) 

Next, we show a method of correction of parallel PSF, it is the same as correcting 
absolute value of the real part. 



Form eq. fl27j) . eq. fl28|) . eq. fl36|) and eq. fl4T|) . we obtain the correction of parallel PSF as 
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follows, 



QsZ^ilfd) 



+Qsq52 



+QsQS2* 



Z'o - ^ (4^0^ - 56* Z',) + ^ {Z', - 26* Z', + 6*'Zt) 



2a4 
^A-^Z!) + ^^{Zt-26Z^,+6^Z^,) 



(///, 5) 



(If S) 



(51) 



Zo^(J^5) ^ g^^oi^/,'^) 



+Q5g52 



^ (2^2* - 35*^2) + ^ (^2^* - 25*^0^ + 6*^Z^) 



-Q5qs2* 



2(t2 

1 

~2^ 



2a^ 
{2Zl - ?>6Zl) + ^ (Z^^ - 25Zo^ + 6'zr) 



(If S) 
(I//, 6). (52) 



Because these equations have only real parts, we obtain 



H'^{I\Zt,,6) = 6 ^ HiiI//,Z^,6) 



+(lsl 



K - ^ ((2 + ^5')H! - 56HI) + ^ ((1 + 56^) H^ - 25(3 + 6-)Zt + (l + 6^) Ht) 
= Hl{I„,Zl6) + Pi{I„,6)qsl 



{I//,Zl6) 



(53) 



When we apply the above equation to a star, we define a distorted image of star /*. 
This image is obtained by using delta function for eq. flT^ ) Therefore I* [6) = P{6) and 
Hl{I* , Zq, 6) = 6. Thus qs2 is obtained as 



qsl 



6-Hi{I*,Zl6) 



(54) 



Therefore, the absolute value of the complex distortion is obtained as 
6 = Hl{I//,Zl6) + ^^^{6-Hl{IlZl6)) 

- Hlih Zl 5) + pp^p^l {5 - Hiil^, Zl 6)) 



(55) 



or 



m{hzi6) 



Pi^il^'^S) rT2 



pE^Iobs,^Sy^2\p/ 



Hl{IhZl6) 



1- 



P^ {I ">">'■, 5) 



(56) 



This equation is the same as the isotropic PSF correction discussed in part 1, because the 
isotropic PSF does not have the orthogonal part. 
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Finally, we obtain the complex distortion in the two dimensional form as follows. 






S = Hl{I„,Zl8)+ X^^^: >^ {S-Hl{I^,Zl6)) 



pEfjobs S\ 

Re [Ul{r'% Zl 5)e-<^^] e^^^ + ^^^^^^.^J ^^ [5 - Re [nl{r''\ Zl 5)6"^^^] e'f^\^l) 



where (ps satisfies as follows 



pEtjobs s\ 

Im [Ul{r'\ zl 5)e-^*^] = p|(^,,J^y m [Ul{r'- , Zl 5)6"^^^] . (58) 



3.5. star selection 

In real analysis, because the ellipticity used in the weight function is obtained by PSF 
correction, we must use some technique to determine it, for example by iteration. In each 
step of iteration, we must determine the centroid and measure the moments of galaxies 
and stars until the result converges, and thus E-HOLICs method need much longer time to 
measure the shear than KSB method. If we use many stars for PSF measurement, the time 
required for the analysis will increase in proportion to the number of stars. This is not an 
essential problem, nut may became a practical difficulty to apply E-HOLICs for wide fields 
surveys. 



4. STEP2 simulation test 

In this section, we perform tests of E-HOLICs method by using STEP2 simulation 
data(Massey et al 2006). This data set contains 6 pattern PSFs, and 64 types of distortion 
fields and the corresponding distortion fields obtained from the original shapes in each PSF 
set by rotating by 90degree. For evaluation of the shear, 1st order polynomial is used to 
estimate the difference between estimated shear and input shear such as 

estimated _ input ^ ^^™P«* _|_ ^ (59) 

where m means over/under-estimation, c means in/over-sufficiency of anisotropic PSF cor- 
rection. More detail information of STEP2 can be seen Massey et al 2006. 
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4.1. star selection 



In real analysis, PSF varies across the field of view, therefore we must use many stars 
for determining PSF distribution. However, because STEP2 simulation data have the same 
PSF in the each fields, we can correct PSF from only one star. To avoid the error from PSF 
measurement, we use only one star which has maximum S/N in each fields. 



4.2. Tests of S/N and size dependence 

A result of STEP2 test in part 1 shows that E-HOLICs method can avoid systematic 
error depended on intrinsic complex distortion | ^«"*''*"^*'^ | _ Thus we restrict ourselves to the 
dependencies of signal to noise ratio "S/N" and size(half light radius "rh") in the shear 
measurement in this paper. 

Figure O shows the results of STEP 2 test by using objects which are S/N>30, 
rh>rhs, 2.5, 3.0 and 3.5 in PSF-A, PSF-B, PSF-C, PSF-D, PSF-E and PSF-F(here and 
after "ABCDEF set") Figure Sis same with Figure [3, but using only PSF-A, PSF-B, PSF- 
C and PSF-F(here and after "ABCF set") where "rhs" is maximum size of half light radius 
of star in each PSF set. These results show that there is a slight overestimation when small 
objects are used, but the overestimation is less than 2a. Figure |5] and Figure |6] are the results 
of STEP 2 test by using objects which are S/N>30, 20, 10, 5 and 3 and rh>3.0 in ABCDEF 
set and ABCF set, respectively. Figure [7] and Figure [H] are the results of STEP 2 test by 
using objects which are S/N>30, 20, 10, 5 and 3 and rh>rhs in ABCDEF set and ABCF 
set, respectively. These results show that there is an underestimation when low S/N objects 
are used, particularly the underestimation is obvious when objects with S/N < 10 are used. 

The details of some results of above selections are listed in Table 1, 
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PSF set 


S/N 


rh 


n(arcmin ^) 


m 


c 


ABCDEF 


30 


3.0(pix) 


2.32878 


0.000634629 


0.0000876483 


ABCDEF 


30 


rhs 


3.85179 


0.00265413 


0.000205060 


ABCDEF 


5 


3.0(pix) 


8.38250 


-0.0212137 


-0.0000551444 


ABCDEF 


5 


rhs 


21.9318 


-0.0255264 


0.000403256 


ABCF 


30 


3.0(pix) 


2.08947 


-0.0045829 


0.000154808 


ABCF 


30 


rhs 


3.68113 


-0.00167376 


0.000225619 


ABCF 


5 


3.0(pix) 


7.84198 


-0.0223220 


-0.0000304242 


ABCF 


5 


rhs 


21.7678 


-0.0244922 


0.000173477 



Table 1. The results of STEP 2 test. In this table, S/N and rh mean the lower limit of S/N 
and rh of used objects, and n means number density of used objects. 

Next, Figjni Fig{TOl Fig{TT] and Fig|T2] show the results in the case of each PSF set, 
respectively. Here the bars with square sign is the result of error for the direction 1 component 
of the shear and the bars with cross sign is for the direction 2 component of the shear. The 
red, green, purple, pink, blue and black colors correspond to PSF-A, PSF-B, PSF-C, PSF-D, 
PSF-E and PSF-F, respectively. We can see PSF dependence from these figures, however, 
there is not only PSF dependence, but also other dependences, and the results are affected by 
combination of them. To clearly see the effect by SN and size, we choose above parameters. 
These results show clearly the underestimation depending on S/N in all PSF set and the 
underestimation in each components in the shear can be also seen by comparing results of 
each PSF sets. 

We can see there is an error in m in psf-B in contrast to the case of psf-A. Data set 
of psf-A and psf-B have same PSF, but background objects have different profiles in these 
sets. The profiles of background objects in psf-B are pure exponential. We guess that the 
difference in the error in m in these data set is due to the difference of radial profile because 
pure exponential is not realistic profile. We can also see there is much larger c error for the 
direction 2 in PSF-D case and for the direction 1 in PSF-E case than other cases, but large 
ellipticities of input PSF are in the direction 1 in PSF-D and in the direction 2 in PSF-E, 
respectively. Therefore these results look rather strange, and similar tendency is obtained 
by other studies (Figure 5 of Massey et al 2007), but only "MJ" in the figure obtained good 
results. Thus there is a possibility that these error may not be due to the method, but due 
to the data itself, and more careful study about PSF correction may be necessary to reveal 
the apparent contradiction. 
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4.3. Summary of STEP2 test 

From the above results of STEP 2 test, E-HOLICs method underestimates the lensing 
shear from low S/N objects, and can estimate the lensing shear precisely if we use only high 
S/N objects or use an appropriate statistical weight which is function of S/N. Systematic 
error depended on size is smaller than that of S/N. This means effect of PSF is smaller than 
random count noise. We do not know a fundamental reason of the S/N dependence found 
above. We need more studies for avoiding and/or correcting this systematic error for the 
precise determination of weak lensing shear (for example cosmic shear study). 



5. A370 galaxy cluster analysis 

In this section, we apply E-HOLICs method to the real data. Our data is a massive 
compact galaxy cluster Abell 370 at z = 0.375 where many distorted images and arcs are 
found. Recently Abell 370 has been studied, for example, Richard et al. 2010 where two 
peaks in central region and elliptical mass distribution in outer region are found by strong 
lensing analysis, Umetsu et al. 2011 and Medezinski et al. 2011 analyzed this cluster by 
weak lensing using Subaru telescope where the color information is used to make a clear 
separation of member galaxies from background galaxies. They are interested in the accu- 
rate determination of the radial mass profile of A370. The data they used was observed 
with the wide-field camera Suprime-Cam (Miyazaki et al. 2002) at the prime focus of the 
8.3m Subaru telescope, and is publicly available from the Subaru archive, SM0KA5. Subaru 
reduction software (SDFRED) developed by Yagi et al. 2002 was used for flat fielding, in- 
strumental distortion correction, differential re-fraction, sky subtraction and stacking. More 
detail information of the data can be seen in Medezinski et al. 2011. We also use the same 
catalogue of background galaxies with Umetsu et al. in this analysis. Medezinski et al. 2010 
and Medezinski et al. 2011 also analyzed A370 with weak lensing method. 



We used IMCAT(http://www.ifa.hawaii.edu/~ kaiser/imcat) and perl language scripts 



(K. Umetsu 2006, private communication) for detecting objects. Stars for PSF correction 
are selected with parameters of 1.60 < rh < 1.95(pix), 19.0 < MAG < 21.0 and 4 < SN, 
and fig.f lTS]) shows their ellipticities measured by KSB quadrupole moments. Typically, they 
have ellipticities around 0.02 along y-axis. We have used the same sample of background 
objects selected by using color information provided by Umetsu(Umetsu et al. 2011), and the 
parameters of rh > 1.95(pix) ^ 0.4(arcsec) = maximum rh of stars, 25.5 > MAG > 20.0 and 
SN > 5. We used Fourier transformation (see part 1) to transform the shear distribution to the 
mass distribution, where the shear distribution is smeared by Gaussian weight with O.larcmin 
scale. FiglH] is the results of 2 dimensional mass reconstruction of A370 by E-HOLICs 



-16- 



method. For comparison we show also 2 dimensional mass distribution obtained by KSB 
method in FigJTSl E-HOLICs method uses 803 objects which corresponds to 7.872(arcmin~^). 
In FigHU minimum contour is k = 0.1, the width between nearby contours is Ak = 0.1, 
and the peak value of k is 1.147 and S/N is 9.034 where we use RMS of B-mode as noise. 
KSB method uses 776 objects which correspond to 7.067(arcmin~^). In FigJT5| contours are 
similarly withdrawn as FiglT^ but the peak value of k, is 1.156 and S/N is 9.788. 

We can see double peaks and elliptical mass distribution in E-HOLICs method which is 
consistent with the result of strong lensing(Richard 2010). On the other hand, KSB method 
cannot detect these peaks. We find that the result of KSB method has no background objects 
just behind A370, but E-HOLICs method has. This will be the reason KSB method can not 
resolve double peaks of mass distribution. 



6. Summary and discussion 

It is widely recognized now that there are many sources of systematic errors in the 
evaluation of weak lensing shear and their improvement is urgently required for the precise 
measurement of cosmic shear. We have developed a new method of shear estimation, based 
on the KSB method and called "E-HOLICs method", by introducing the elliptical weight 
function to define multipole moments of galaxy light distribution in previous paper (part 1). 
The use of elliptical weight function is expected to avoid the systematic error coming from an 
expansion of the weight function, which was usually done in some of the previous approaches 
to shape measurement, including the KSB method. Bernstein & Jarvis 2002 and Melchior 
2011 are also developing schemes which use an elliptical weight for the shape measurement. 

Following part 1 we developed the correction scheme of PSF with anisotropy in E- 
HOLICs method. Generally, the direction of the elliptical PSF is not same with the direc- 
tion of distortion, and thus we must treat these two direction in PSF correction. By dividing 
anisotropic part of PSF into parallel and orthogonal part against the direction of the distor- 
tion, we succeeded to reduce PSF correction into two steps of simple one dimensional form. 
Then we test E-HOLICs using STEP 2 data simulation. In particular we tested the depen- 
dence on the S/N and size of background galaxies. We have found that E-HOLICs gives 
underestimation for low S/N sources. Although the precision is of the same level of other 
methods, this causes a serious problem in the accurate shear measurement. We have found 
that at least a part of reason for the underestimation is the errors in the centroid. Although 
the error in the position of centroid distributes randomly, the area in the image where the 
measured ellipticity from the centroid within it is smaller than the correct value is larger 
than the other area in the image. Thus the random error in the centroid has a tendency to 



-17- 



cause the underestimation of the measured shear. We are investigating this effect and the 
result will be shown in the forthcoming publication. 

In this paper we have ignored higher-order moments in PSF which might be important in 
some cases. We aim to study the higher-order PSF correction also in the future publication. 

As we pointed out, E-HOLICs method requires longer time for measuring moments if 
we use many stars for PSF correction than other methods such as KSB method. This might 
cause a practical difficulty in any wide field surveys. We hope that there will be a way to 
avoid this problem by using powerful computer specially programed for this purpose. 

As an application of E-HOLICs methods, we analyzed a massive and compact galaxy 
cluster Abell 370 using by Subaru/S-Cam data. We found that E-HOLICs can detect two 
mass peaks which is consistent with the strong lensing analysis. The KSB method cannot 
detect two peaks possibly because available number of background objects is reduced in 
compared with E-HOLICs. 

Although the E-HOLICs method developed in this paper has a potential to accurately 
measure the shear, it has still some shortcoming as described above and definitely more stud- 
ies are necessary to apply E-HOLICs to the planned large-scale cosmic shear observations. 

We thank K. Umetsu for a useful discussion and for providing his scripts, N. Okabe for 
useful discussion. This work is supported in part by the GCOE Program "Weaving Science 
Web beyond Particle-Matter Hierarchy" at Tohoku University and by a Grant-in-Aid for 
Scientific Research from JSPS(nos. 18072001, 20540245 for TF), as well as by Core-to-Core 
Program "International Research Network for Dark Energy." 
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A. List of newly defined functions 



Function 


Definition from 


Ellipricity 


Name 


i^CP) 




circular 


zero image 


P'°0) 


= Jd^(i)l^{4))P''°0-4)) 


circular 


smeared zero image 


iHo) 


= i'CP) 


having a ellipticity of S 


lensed image 


Jobsf^Q^ 


= ^d^i^{<p)p{e-<i^) 




observed image 


m 


= j (mi\<p)p,^e - <p) 


having a direction of distortion 


parallel image 


He) 


= Jd'<i)i^i<l>)p{e-(t>) 


having a ellipticity of S 


distortional image 


Pie) 


= J d'<PP0)q,i9 - <p) 




real PSF 


P^'\J3) 


= p{e) 


circular 


isotropic zero PSF 


m 


= J d^P//{<P)qs{e - <p) 


having a ellipticity of S 


distortional PSF 


p/m 


= jd^(t>P''"{<P)q/ke-<P) 


having a direction of distortion 


parallel PSF 


qM 


P{9) = f d'<pp/^<p)q,ie - <p) 


orthogonal to distortion 


orthogonal part of P{9) 


QsiO) 


P{9) = fd'^P//{<P)qs{9-<p) 


having a direction of distortion 


distortional part of P{9) 
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Fig. 1. — Relations between zero image, source and image. Black arrow shows calculation 
in E-HOLICs methods 
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Fig. 2. — The relations in PSF correction. Black arrows show calculations in E-HOLICs 
methods 
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Fig. 3. — Result of testing ABCDEF set by using objects which are S/N>30 and rh>rhs, 
2.5, 3.0 and 3.5 
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Fig. 4. — Result of testing ABCF set by using objects which are S/N>30 and rh>rhs, 2.5, 
3.0 and 3.5 and 
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Fig. 5. — Result of testing ABCDEF set by using objects which are S/N>30, 20, 10, 5 and 
3 and rh>3.0pix 
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Fig. 6. — Result of testing ABCF set by using objects which are S/N>30, 20, 10, 5 and 3 
and rh>3.0pix 
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Fig. 7. — Result of testing ABCDEF set by using objects which are S/N>30, 20, 10, 5 and 
3 and rh>rhs 
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Fig. 8. — Result of testing ABCF set by using objects which are S/N>30, 20, 10, 5 and 3 
and rh>rhs 
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Fig. 9. — Result of each PSF set and each direction by using objects which are S/N>30 
and rh>3.0pix 
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Fig. 10. — Result of each PSF set and each direction component by using objects which are 
S/N>5 and rh>3.0pix 
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Fig. 11. — Result of each PSF set and each direction component by using objects which are 
S/N>30 and rh>rhs 
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Fig. 12. — Result of each PSF set and each direction by using objects which are S/N>5 
and rh>rhs 
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Fig. 13. — Elhpticity of stars which are used for measuring PSF. These eUipticity are 
elongated by 10000 times, so roughly they have elhpticity 0.02 along y-axis. 
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Fig. 14. — A370z analyzed by E-HOLICs method. Contours mean A/t = 0.1 and minimum 
contour is 0.1. 
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Fig. 15. — A370z analyzed by KSB method. Contours mean Ak 
contour is 0.1. 
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